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We revisit a simple model that combines minimal gauge mediation and the next-to-minimal 
supersymmetric standard model. We show that one can obtain a 125 GeV Standard Model-like 
Higgs boson with stops as light as 1.1 TeV, thanks to the mixing of the Higgs with a singlet state 
at 0(90 — 100) GeV. Sparticle searches at the LHC may come with additional 6—jets or tans and 
may involve displaced vertices. The sparticle production cross-section at the 13 TeV LHC can be 
0(10 — 100) fb, leading to great prospects for discovery in the early phase of LHC Run H. 


I. INTRODUCTION 

The discovery of a Higgs-like scalar particle with mass 
close to 125 GeV [T] has considerable impact on super- 
symmetric (SUSY) model building. In its simplest re¬ 
alization, the Minimal Supersymmetric Standard Model 
(MSSM), the tree-level Higgs mass is bounded from 
above by the Z-boson mass, which implies that large 
radiative corrections of the order of the tree-level mass 
are needed [2]. This motivates extensions of the mini¬ 
mal model with new tree-level contributions to the Higgs 
mass. A possible source of enhancement of the tree- 
level Higgs mass is mixing with an additional neutral 
state that is lighter than the SM-like Higgs. This situa¬ 
tion can be realized in the Next-to-MSSM (NMSSM) [3]. 
Analyses of the generic NMSSM parameter space have al¬ 
ready demonstrated that this possibility is viable HE]- 
Here instead we want to study this scenario in a simple 
and predictive framework of SUSY breaking, gauge me¬ 
diation [5], which solves elegantly SUSY CP and flavor 
problems. 

Indeed the combination of the NMSSM and gauge me¬ 
diation is particularly motivated, as the NMSSM pro¬ 
vides a simple solution to the notorious problem [7] 

of gauge mediation. Yet it is very difficult to realize this 
scenario with minimal gauge mediation (MGM), as the 
NMSSM soft terms are too small [5]. These problems 
can however be cured by adding direct couplings of the 
singlet to messengers, at the cost of a single new pa¬ 
rameter (for a different approach see e.g. Ref. |^). A 
viable model of this kind has been proposed by Delgado, 
Giudice and Slavich (DGS) in Ref. [10] (extensions with 
Higgs-messenger couplings were studied in Ref. my 
However, the authors of Ref. m concluded that in this 
model sparticles cannot be lighter than in MGM. 

In this article we re-analyze the DGS model and iden¬ 
tify new viable regions in the parameter space where 
singlet-Higgs mixing is small enough to pass experimen¬ 
tal constraints, but large enough to give substantial con¬ 
tributions to the tree-level Higgs mass. This model can 


therefore rely on smaller contributions from stop loops, 
thus reducing the overall scale of sparticle masses. In¬ 
terestingly, squarks and gluinos can be light enough to 
be discovered in the early stage of the LHC run H, in 
contrast to MGM, where a 125 GeV Higgs mass requires 
colored sparticles beyond the reach of the LHC (even 
for very high luminosity) |12j . Moreover, we find that 
the light singlet-like scalar can easily explain the 2a ex¬ 
cess around 98 GeV observed in the LEP Higgs searches 
[laiiii. The realization of this scenario, with maxi¬ 
mal contribution to the tree-level Higgs mass from mix¬ 
ing, hxes almost all of the model parameters. A single 
parameter remains free and controls the details of the 
phenomenology. The lightest supersymmetric particle 
(LSP) is the gravitino and the next-to-LSP (NLSP) is 
the singlino, a setup that leads to new signatures at col¬ 
lider experiments. The underlying model might therefore 
serve as a representative for a whole class of signatures 
that motivate suitable SUSY search strategies. 

II. THE DGS MODEL 

The field content of the DGS model (see Ref. [10] 
for details) consists of the NMSSM fields (the MSSM 
fields plus a gauge singlet S), in addition to two copies 
of messengers in 5 - 1-5 of SU(5), denoted by 
i = 1,2 with SU(2) doublet and SU(3) triplet compo¬ 
nents <l>f, , i = 1,2. Supersymmetry break¬ 

ing is parametrized by a non-dynamical background field 
X = M -\- FO^. Apart from the Yukawa interactions, the 
superpotential is given by the NMSSM part, the spurion- 
messenger couplings and the singlet-messenger couplings, 
W = IUnmssm + bPcM + where 

kkNMSSM = XSHuHd -\- -S^ , (1) 

VEgm = ^ X , (2) 

Wd^GS = S {^D^1^2 ■ 


( 3 ) 
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The new couplings in Eq. ([^ are assumed to unify at 
the GUT scale ^d(Mgut) = Ct(-Wgut) = ^. A similar 
assumption can be made for but these parameters 

are largely irrelevant for the spectrum. 

Z 3 invariant NMSSM models such as this one have 
a potential cosmological problem with domain walls, 
which are predicted to appear during the phase transition 
associated with electroweak symmetry breaking. This 
can be solved by the introduction of higher dimensional 
Za-violating operators. Dangerous tadpoles that may 
threaten successful electroweak symmetry breaking may 
be avoided by the imposition of a discrete 7?.—symmetry 
on the higher dimension operators |15|. 

Through the superpotential in Eq. Q the messengers 
feel SUSY breaking at tree-level and communicate it to 
the NMSSM fields via gauge interactions and the direct 
couplings in Eq. ([^. The contribution from gauge inter¬ 
actions is given by the usual expressions of MGM for one- 
loop gaugino masses Mi and two-loop sfermion masses 
mj at the messenger scale M 
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Mi = 2gfm , = 4 ^ C^{f) , (4) 

i=l 


where m = l/(167r^)E/M and Ci{f) is the quadratic 
Casimir of the representation of the field / under SU(3) x 
SU(2) X U(l). The contributions from direct singlet- 
messenger couplings generate one-loop A-terms for the 
NMSSM couplings 

Ax = ^=-m {2^1 + ^1) , (5) 


and two-loop contributions to soft masses for the singlet 
and the Higgs fields 

rh% = ifi^ -I- 15^1- -I- 

[4^2 (2^1^ ^3^2)] 


(2^1) -f . 


( 6 ) 

( 7 ) 


There is also a one-loop contribution to the singlet soft 
mass nni that is relevant only for very low messenger 
scales and has negligible impact on the spectrum. The 
model is thus determined by five parameters: m, M, A, k 
and where one parameter (following DGS we choose n) 
can be eliminated by requiring correct electroweak sym¬ 
metry breaking (EWSB). 


III. LOW-ENERGY SPECTRUM 

There are several regions in the parameter space that 
lead to a viable spectrum and are compatible with per¬ 
turbative couplings up to the GUT scale. They can be 
distinguished by the size of the relative contributions to 


the SM-like Higgs mass, which are given schematically 
by 

= Aff cos^ 2j3 + sin^ 2/? -f + '^'h^ioop j (8) 

where m\ is the contribution from mixing with the 
other two CP-even states in the full Higgs mass matrix. 

Since a larger tree-level mass implies a lighter SUSY 
spectrum, we concentrate here on the region in param¬ 
eter space where the effective tree-level contribution to 
the Higgs mass is maximized. It turns out that one can 
reach up to ml — loop ~ (99GeV)^, provided that 
tan P is large and the contribution from mixing is sizable 
and positive. This requires the singlet state to be lighter 
than the SM-like Higgs, which is not excluded by LEP 
if the mixing angle with the SM-Higgs is small enough. 
Note that the LHC constrains this scenario only through 
measurements of the SM-like Higgs couplings that are 
suppressed by the mixing. 

A focal feature of this scenario is its high level of pre- 
dictivity, as three out of four free parameters of the model 
are determined by the Higgs sector. Maximizing the tree- 
level contribution to the SM-like Higgs mass 771^,2 fixes the 
singlet-like Higgs mass and the singlet-Higgs mixing 
angle 9 (mixing with the heavy Higgs doublet is negli¬ 
gible) along the lines of Ref. [S], giving approximately 
m/ij « 94 GeV and cos 9 ~ 0.88. This in turn determines 
the model parameters A and while the overall scale of 
soft terms m is fixed by the required size of the residual 
loop contribution to mh 2 ■ Departing from the maximal 
tree-level contribution leads to slightly different mhi and 
9, and of course larger m. 

In practice one can map the model parameters 
to the low-energy spectrum only numerically. For 
this analysis we have used a modified version of 
NMSSMTools [16]. Independent checks using a modified 
version of S0FTSUSY3.4.1 m produced Higgs and SUSY 
spectra that agreed at the percent level. 

Before discussing our results, we provide some rough 
analytic results that can be obtained neglecting renor¬ 
malization group (RG) effects and expanding the 
NMSSM vacuum conditions in the limit of large singlet 
vacuum expectation value. In this way one obtains ap¬ 
proximate relations between A and the physical Higgs 
parameters 


ruhi 

A-/2g3rh 




4vrh 


i sin 20. (9) 


For TeV-scale superpartners (and the above values for 
m/ij and 9 that maximize the Higgs tree-level contribu¬ 
tion) one finds ^ ^ A ^ 10“^. 

The smallness of ^ and A has important consequences 
for the low-energy spectrum. Small ^ implies small values 
of |Aa| and \A^\ and imposing proper EWSB yields n 
A and the prediction of large tan/3 ~ A/k. In turn, the 
smallness of and k results in a very light singlet-like 
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pseudoscalar oi of mass 


rua 


3253 


-mhi . 


( 10 ) 


For ^ ^ 10 Eq. (10) predicts to be smaller than 


m/jj by a factor of a few. We find numerically that the 
light pseudoscalar mass varies between 20 and 40 GeV, 
for gluino masses below 2.5 TeV. 

For the above range of parameters the mass of the 
singlino can be obtained from the following approximate 
sum rule m- 


’ (11) 

which implies that the singlino mass is about 100 
GeV. Since in MGM the LSP is the gravitino and the 
typical scale of the NLSP is the bino mass Mi « 
420 GeV (m/TeV), it is clear that here the singlino 
strongly dominates the composition of the NLSP. This 
is a distinguishing feature of this model. 

This is closely connected to the main virtue of this 
scenario, the large contribution to the tree-level Higgs 
mass from singlet-Higgs mixing. This requires smaller 
radiative corrections from stop loops, and in turn much 
lighter sparticle masses than in MGM. Through these 
corrections the observed Higgs mass of 125 GeV essen¬ 
tially fixes the overall scale of the sparticle spectrum to, 
up to an estimated theoretical uncertainty of 3 GeV in 
the prediction of rrih^- We find that mh 2 = 125 GeV is 
compatible with a gluino mass of 2.1 TeV (1.4 TeV if 
the theoretical uncertainty of 3 GeV for the Higgs mass 
is taken into account). Squarks of the first two gener¬ 
ations have approximately the same mass as the gluino 
and should be within the reach of the LHG Run H. Stop 
masses can be as light as 1.7 (1.1) TeV, which should be 
compared with the lower bound on stop masses of about 
8 (3) TeV in MGM [T^]. In Fig. [^we show the values of 
Higgs, gluino and stop masses for several model points 
separated by the singlet-Higgs mixing angle 6 . Note 
that cos^ 9 is roughly of the size of effective Higgs sig¬ 
nal strengths Ri = {a x BR)i/(CT x BR)®'^, which are 
substantially reduced in this scenario. Nevertheless all 
shown points are compatible with LEP and LHG con¬ 
straints on the Higgs sector. 

Having fixed (^, A, to) by the set of physical Higgs pa¬ 
rameters (TOft,j, TOh 2 , 6 *), the only free parameter left is 
the messenger scale M. This parameter controls the low- 
energy spectrum in several ways. First of all, increasing 
M leads to larger values of At at the electroweak (EW) 
scale, which (as in MGM) is purely radiatively generated 
and therefore grows with the length of the RG running. 
In turn, this enhances the stop-mixing contribution to 
the Higgs mass, and therefore larger M leads to lighter 
stops and hence smaller to for the same value of mh 2 - 
Also, the value of M essentially determines the nature 
of the next-to-NLSP (NNLSP). For small M < 10® GeV 
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FIG. 1: (Upper panel) SM-like Higgs mass vs. gluino mass 
and (Lower panel) singlet-like Higgs mass vs. lightest stop 
mass. The various model points are distinguished by the 
Higgs-singlet mixing angle 6 , which decreases from top to 
bottom as specihed in the upper panel. For the same SM- 
like Higgs mass a larger mixing angle allows for much lighter 
gluinos. The lightest stop masses are obtained for a singlet- 
like Higgs around 94 GeV. 


the (mostly right-handed) stau is the NNLSP (with se- 
lectron and smuon being co-NNLSP), because the soft 
mass TO^ is smaller than Mi at the messenger scale. For 
M > 10® GeV (requiring gluino masses below 2.5 TeV) 
the RG effects are strong enough to raise to^ above Mi 
and the bino-like neutralino becomes the NNLSP. In the 
transition region 10® GeV < M < 10® GeV the NNLSP 
can be either stau or bino, depending on the other pa¬ 
rameters. The messenger scale M controls the gravitino 
mass according to: 


TO 3/2 = 38 eV 



M \ 
10® GeV ) ■ 


( 12 ) 


The simultaneous presence of gravitino LSP and singlino 
NLSP leads to a novel phenomenology quite different 
both from MGM models and from typical NMSSM sce¬ 
narios. 
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IV. LHC PHENOMENOLOGY 


In our scenario decay chains of every supersymmetric 
particle produced at the LHC end up in a singlino-like 
neutralino Ni. Since the singlet couples very weakly, 
these decays always proceed through the NNLSP or co- 
NNLSP. The singlino subsequently decays to the grav- 
itino and the singlet-like pseudoscalar ai, which in turn 
predominantly decays to 6-quarks: 

Ni^aiG^bbG. (13) 


The decay length of the neutralino (in its rest frame) is 
approximately given by 


‘Ni 


2.5 cm 


/100 GevV 


M 

106 GeV 



(14) 

Since M cannot be much below 10® GeV, it is clear from 
the above formula that the singlino NLSP (with mass 
about 100 GeV) always travels macroscopic distance be¬ 
fore it decays. For large M the singlino decays well out¬ 
side the detector so it is stable from the collider point 
of view. However, for M ~ 10® — lO’’ GeV the singlino 
may decay in the detector after traveling some distance 
from the interaction point leading to a displaced vertex. 
Since the value of M also decides about the nature of 
the NNLSP, it can be used to define three regions with 
distinct LHG phenomenology, which we briefly discuss in 
the remainder of this letter. A more detailed analysis of 
LHG phenomenology and discovery prospects will be the 
subject of a future publication. 

In Table U we collect several characteristic benchmark 
points. Points PI and P4 represent the lightest SUSY 
spectra we have found, for very low and very large mes¬ 
senger scales, respectively. Since the Higgs mass errors 
are pushed to the limits, we consider these points merely 
as limiting cases, although not necessarily unrealistic. 
Note in particular that P4 is not obviously ruled out by 
standard SUSY searches for jets -I- missing Et, since the 
additional decay of the would-be-LSP bino to singlino re¬ 
duces efficiency compared to the GMSSM [T51 . The 

other points are representatives for the three characteris¬ 
tic regions discussed below, and P3 is in addition cho¬ 
sen to fit the LEP excess. Note that all points have 
quite large singlet-Higgs mixing, leading to reduced effec¬ 
tive Higgs couplings. Points with smaller mixing and/or 
larger Higgs masses can be obtained by increasing the 
overall SUSY scale rh. 

In all regions sparticles can be very light, so that huge 
parts of the parameter space are in the reach of LHG Run 
H. As can be seen from Table |l] the total strong produc¬ 
tion cross-section (dominated by qq and qg) is 0(10—100) 
fb, as computed with PROSPIND [^. LHC Run H is ex¬ 
pected to deliver 0(10) fb“^ of integrated luminosity in 
2015, which results in 0(100 — 1000) potentially discov¬ 
erable events. The total EW production cross-section at 
the 13 TeV LHC (computed with Pythia 8.2 [3T]) is 



PI 

P2 

P3 

P4 

P5 

m 

7.5 • 10 

^ 8.7-10^ 

9.3- 10^ 

5.9- 10^ 

9.3- 10^ 

M 

1.4- 10 

® 2.8-10® 

3.3- 10'^ 

8.3 ■ 10^"^ 

3.4- lO^'^ 

A 

1.0- 10“ 

2 9.3- 10“^ 

6.7- 10"® 

9.2 ■ 10“® 

6.9- 10“® 

c 

1.2 • 10~ 

2 1.1 • 10~^ 

1.3- 10"^ 

3.2- 10“^ 

2.0- 10“^ 

tand 

25 

28 

24 

26 

21 

mhi 

92 

93 

98 

94 

94 

mh2 

122.1 

123.4 

122.9 

122.1 

125.0 

rriai 

26 

26 

28 

40 

32 


101 

102 

106 

104 

104 


322 

377 

400 

251 

379 


303 

358 

406 

449 

676 


284 

333 

376 

432 

637 

m-g 

1.73 

1.98 

2.09 

1.37 

2.06 

rriuj^ 

1.79 

2.06 

2.15 

1.36 

2.07 

■mi, 

1.64 

1.87 

1.90 

1.06 

1.63 


6.4- 10“ 

2 0.34 

48 

1.9 ■ 10^® 

6.0 ■ 10^® 

_13TeV 

^qq 

9.35 

2.99 

1.98 

59.7 

2.63 

_13TeV 

11.9 

3.30 

2.01 

91.1 

2.48 

13TeV 

^strong 

25.2 

7.28 

4.58 

190 

5.95 

8TeV 

strong 

0.51 

0.07 

0.03 

10.1 

0.05 

13TeV 

^EW 

27 

12 

7.5 

6.7 

5.6 

^8TeV 

^EW 

5.5 

2.1 

1.2 

1.3 

0.7 


TABLE I: List of benchmark points. All masses are in GeV 
except colored sparticle masses in TeV, the neutralino decay 
length in m and cross-sections in fb. All points have 
reduced effective Higgs couplings, with Higgs signal strenghts 
about 0.75, as a result of a Higgs-singlet mixing angle with 
cos 6 « 0.88. 


typically comparable to the strong one but is distributed 
among many different channels with rather small indi¬ 
vidual cross-sections of order 0(1 — 10) fb. The most fre¬ 
quent EW production channel is xt Xs (which are wino- 
like states decaying dominantly to staus) with the cross- 
section of about one fifth of the total EW cross-section 


A. Low-M Region: M < 10^ GeV 


In this region, represented by benchmarks PI and 
P2 in Table |Tj the lightest stau is the NNLSP (with 
smuon/selectron co-NNLSPs) and therefore the singlino 
is produced in association with either tau or leptons. 
Since the splitting between sleptons and the singlino is 
around 200 GeV or more, one expects high-py taus or 
leptons in the final state, which presumably can be used 
to reduce QGD backgrounds considerably. In this region 
the singlino decays (via light pseudoscalar) to 66 still in¬ 
side the detector. However, identifying these displaced 
6 -jets might be challenging since they are expected to be 
very soft due to the small pseudoscalar mass. We note 
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that the low-M region is constrained (as are all GMSB 
models) by the matter power spectrum as inferred from 
the Lyman—a forest data and WMAP [^, which dis¬ 
favours a gravitino mass between 16 eV and where 
O(keV). In fact, from our scan we do not find any 
gravitinos with mass less than 16 eV. The other bound 
implies that 


cTfi, ^ 17 m 


100 GeV 


M 


Ni 


'"' 3/2 


(15) 


from Eqs. (12|,(|T4|. We remind the reader that is 
close to lOOGevbecause of the sum rule Eq. (11). PI 


and P2 violate this bound, whereas P3 may or may not, 
depending on the precise value of rriy'^- However, en¬ 
tropy production after gravitino decoupling may provide 
a cosmological evasion of the bound for any point 


B. Medium-M Region: lO'^ GeV < M < 10® GeV 


V. SUMMARY AND CONCLUSIONS 


In this letter we have re-analyzed the DGS model, a 
simple and predictive framework for combining MGM 
and the NMSSM. We have found new regions in the pa¬ 
rameter space with a singlet at 0(90 — 100) GeV and 
singlet-Higgs mixing giving substantial contributions to 
the tree-level Higgs mass. While these regions are com¬ 
patible with Higgs precision data, we find that colored 
sparticles can be close to their direct search limits, in 
sharp contrast to MGM. 

SUSY decays have more visible particles and can lead 
to less missing Et as compared to MSSM predictions. 
The phenomenology is controlled by a single parameter, 
which determines whether SUSY decays chains lead to 
additional 6—jets or taus and involve displaced vertices. 
As the production cross-sections of colored sparticles are 
0(10 — 100) fb, significant parts of parameter space are 
discoverable in the first year of LHG Run H. 


In this region, represented by benchmark P3, the 
singlino LSP is long-lived. Stau is still NNLSP, but 
smuon/selectron are no longer co-NNLSPs because they 
are heavier than the bino-like neutralino. In consequence, 
a vast majority of gluino and squark decay chains ends 
in stau NNLSP decaying to tau and quasi-stable singlino 
NLSP, with two high-pr taus in each event. 

C. Large-M Region: M > 10® GeV 

Eor large messenger scales, represented by benchmarks 
P4 and P5, the NNLSP is bino-like. Therefore the (quasi¬ 
stable) singlino is typically produced in association with 
the 125 GeV Higgs, BR(iV 2 ^ iVi/ 12 ) ~ 70 - 75%, or 
the singlet-like Higgs, BR(iV 2 — Nihi) ^ 25 — 30%. 
Both Higgs states decay dominantly to bb. Using a 6—jet 
tagging efficiency of 70% [IS], one still expects in each 
event at least two (three) identified high-pT 5—jets from 
bino decays with a probability of about 60 (30)%. This 
comes on top of the 6—jets originating from other decays 
in the gluino and/or squark decay chains. Therefore, this 
model can be easily discriminated against MSSM models 
using searches with large numbers of 6—jets. 
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